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We reexamine critically the chiral expansion for the baryon magnetic mo- 
ments including the contributions from loops which involve intermediate octet 
and decuplet baryons. We find that, contrary to some claims, the nonanalytic 

1/2 

loop contributions of orders ms and rus In mg are of the same general size 
because of large coupling factors for the latter, and that the decuplet con- 
f\J I tributions are as large as the octet contributions and must be included in a 

2l ! consistent calculation. There is no clear evidence of the convergence of the 

^Sj ' chiral series. The adequacy of the theory will not be established until dynam- 

ical models are able to calculate the contributions from the counterterms that 

t — ■ 

Qs . largely hide the loop effects in fits to the data. 
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I. INTRODUCTION 

Chiral perturbation theory (ChPT) has been a very useful approach to the theory of 
the low momentum processes involving mesons and baryons, and has been used in various 
attempts to explain the baryon magnetic moments. All the moments in the baryon octet 
except that of the S*^ have been measured. The S'^A transition moment /xsa is also known. 
At the SU(3)-symmetric tree level [ij, ChPT parameterizes these eight measured quantities 
terms of the two parameters /i^i and fip in the effective Lagrangian 

C = -^{fiDTr Ba^^F^^^iQ^B} + fipTi Ba^^F^^iQ^B]) , (1) 

and is able to fit them with an average error of about 0.24 nuclear magnetons /iN- In 
this expression, B is the usual representation for the baryon-octet in flavor-space, Q = 
diag(2, —1, — 1)/3 is the quark matrix, m^ is the nucleon mass, and F'^'^ is the electromag- 
netic field. 

The study of the baryon moments beyond the tree level has a long history which can 
be found in [^,p|,^,p|,P,|7|j8|j9|,p!0|,p!l[] . The approaches and conclusions have varied. Caldi and 
Pagels pointed out that the loop corrections lead to uniquely defined "nonanalytic cor- 
rections" to the moments proportional to the square root and logarithm of the symmetry 
breaking parameter rris of ChPT. In contrast, the analytic corrections involving integer 
powers of rris are cutoff dependent in perturbation theory. The contributions left after the 
divergences are subtracted are not uniquely defined in the absence of further dynamical 
information, and mix with the contributions of higher chiral couplings. Since these contri- 
butions are expected to be small, they have generally been ignored. Caldi and Pagels 
calculated what were supposed to be the leading corrections to the baryon moments, those 
of order ml^'^, and found that they are as large as the tree- level moments and fail to improve 
the fit to experiment. These terms were recalculated by Gasser et al. 0, and the logarith- 
mic terms of type m^Yimg were analyzed by Krause ^ , but without detailed analyses of the 
moment problem. 

The problem has been reexamined in recent years. Jenkins et al. [^ used heavy baryon 
chiral perturbation theory (HBChPT) [|l^] to calculate the one-loop corrections to the octet 
baryon moments allowing both octet and decuplet intermediate states. Their formal analy- 
sis included the nonanalytic corrections of types m^^ and mslnrris and all the counterterms 
allowed to that order. However, they assumed in fitting the moment data that the loga- 
rithmic corrections and the counterterms were small compared to the terms proportional to 
TTig'^, and included only the last in their analysis. The resulting fit to the data was much 
worse than the tree-level fit if they used tree-level values for the axial coupling constants, 



but better if they used the axial coupling constants extracted from one-loop fits [|14| to other 
data, F = 0.4, D = 0.61, and C = —1.2, values about 20% smaller than those derived at tree 
level. This sensitivity leaves the significance of the loop corrections unclear. In contrast. 



Meissner and Steininger [|T0| considered all octet terms (analytic and nonanalytic) which are 



of order g^ or less in ChPT. They did not include the intermediate decuplet states directly, 
but argued that some decuplet contributions are included because of the way they deter- 
mined the couplings. With the counterterms included, they could fit the seven well measured 
octet moments precisely, and predicted a EA transition moment of (1.40 ± 0.01 /xn) to be 



compared with the experimental value, /iea = (1.610 ±0.080) /iN. Their principal conclusion 
was that the chiral expansion appeared to be converging rather well. 

The problem has also been studied using the l/N^ expansion by a number of authors 
@j^jl)i] who showed that the m^^ terms are leading order in the expansion, and obtained a 
number of interesting sum rules for the moments. Dai et al. again kept only the leading 
nonanalytic corrections proportional to my^ and the counterterms in their analysis of the 
moment data, and obtained a good overall fit, but with seven parameters. Finally, Bos et al. 



Ill considered the moments from the point of view of flavor SU(3) breaking in the baryon 
octet using a chiral counting and treatment of meson masses different from that usually used 
and again obtained a successful parametrization using seven parameters to describe eight 
measured moments. 



II. PREDICTIONS OF THE CHIRAL EXPANSION 

Thus far, the numerically successful approaches to the baryon moments have all used 
seven parameters to fit eight experimental data. Two of these parameters are the tree-level 
moment parameters ^d = ^i and fip = ^2. The remainingfive are the coefficients 03, . . . , ay 
of the SU(3) symmetry breaking counterterms defined as HI 

/:,t = -r^F^''{asTTB[[Q,a,,B],M] + a.Ti B{[Q,a,,B],M} 

+a,TicB[{Q, a^^B}, M] + aeTiB{{Q, a^,B},M} (2) 

+a7TrBa^^BTiMQ+ ...) . 

Here A^ = diag(0, 0, 1) is a matrix used to introduce SU(3) breaking introduced by the 
strange quark mass. A factor r/ig has been absorbed in the coefficients of the Lagrangian. 
There are two terms with the same structures as those in the tree level moment operator 
(|I|) hidden in the ellipsis . 

Given the large number of parameters used in the fits, it is natural to ask if the loop 
corrections matter, and if one can show that the uniquely extractable nonanalytic corrections 
actually improve the tree-level fit to the baryon moments. To answer these questions and 
investigate the expected cancellations and the convergence of the chiral expansion, we have 
redone the calculations of Jenkins et al. 0, corrected several errors in the published work, 
and reanalyzed the problem of fitting the measured moments from a different point of 
view. We have followed Jenkins et al. in retaining only the nonanalytic contributions of 
types m^/^ and mglnms in our analysis, but including the effects of the octet-decuplet mass 
splitting throughout. The graphs which produce these contributions are shown in Figs. 1 
and 2, respectively. 



^We follow the notation used in Q except for a rearrangement of the counterterms following Bos 
et al. dnl. 

^There is a factor (—5) missing in a contribution from the graph 2c and a factor 1/2 missing from 
the decuplet magnetic moment operator. Other corrections are noted in the erratum to 0]. 



The octet baryon moments are given in the one-loop approximation by 

where iJ,\ is the tree-level moment, fx\ contains the contributions from loops that involve 
only intermediate spin-1/2 octet baryons, and ^l contains the loop contributions that 
involve the spin-3/2 decuplet baryons, 



nf^ = Oi , (4) 

(1/2) _ Y^ rn^Mx ^(x) 

(X) _^,(X)^_ x,.2i .^2 



1^"' = E ^^S9^A'"' + 

1 



E TT^(7r^-2A,^-^^a.)MilnMi, (5) 

X=-K,K,r) J 

X=TT,K °'^J 

E T^^[(7P-2AP^a.)^(3/2)(Mx,5,/.) +7Pl;3/2)(Mx,5,/.)]. (6) 

=TT.K.n J 



X=TT,K,ri 

The functions F, L, and G are defined as 

^2 f 2VM2 - 52 [;r/2 - arctan (5/VM2 - 52)], M > 5, 

TT F(M, 5, /x) = -5 In— + , (7) 

^ [ -2V(52 _ ]vf2 In [(5 + ^52 _ M2)/M], M<5, 

/Vf2 
L(3/2)(M,<5,/i) = M^ln^ + 27r5F(M,<5,/i), (8) 

/i 

L;3/2)(M, 5, /x) = M' In^ + ^ G(M, S, /i), (9) 

M2 

7rG'(M,5,/i) = -5Hn^ + nM^ 



2(M2 - ,^2)3/2 [^/2 - arctan {6/VM^ - 5'^)], M > 6, 

+ { , ■ (10) 

-2(52 - Af 2)3/2 In [(5 + V52 _ M2)/M], M < 5, 

The coupling coefficients ctj, /3i, /?j, Aj + Aj, 7j + 7j + 7j are identical to those in P where 
Aj + Aj is denoted by Aj and 7i + % + 7i is denoted by 7j. The necessary corrections have 
been included. The counterterm contributions follow from Eq. (^ and are given in Table I. 
Before examining the loop contributions in detail, we note that the seven well-measured 
octet moments can be fitted exactly using the five counterterm couplings 03, . . . , 07 and the 
two tree-level parameters fin = 0,1 and fip = 0,2, even in the presence of the loop corrections-^! 
To see this, we define an 8 x 7 matrix X which has as its ith column the matrix elements of 



^We choose to focus on the seven baryon moments in our analysis because of the high accuracy 
with which they are known relative to /isA- The alternative of making a weighted least-squares 



the corresponding operator in the physical baryon states p, n, E+, E~, S", S~, A and the 
S'^A transition matrix element as in Table I. The most general column vector of magnetic 
moments consistent with this chiral structure is then of the form 

/i = Xa, (11) 

where a is a seven-component column vector of coefficients. The 7x7 matrix X obtained by 
omitting the S°A row in X is invertible, so coefficients Oj can be found which give exact fits to 
either the seven well- measured baryon moments, a = X~^/imeas, or to the difference between 
the measured moments and the loop corrections in this sector, a' = X~^{firacas — Aioop)- 
In the second case, the moments are again reproduced exactly when the loop corrections 
are added in. The only remaining constraint on the theory is that provided by the S°A 
transition moment. 

In the absence of loop corrections, yUsA is determined entirely by the seven coefficients a 
through the product Xa with the full matrix X, or, equivalently, through the Okubo sum 
rule 

4\/3/isA = 6/iA + /is- +/US+ -4yU„-4/iH0- (12) 

The origin of the sum rule can be seen rather simply in the present context. We can extend 
the matrix X to an 8 x 8 matrix X by adding a final column with zero entries without 
affecting the seven baryon moments. That is, fi = Xa with the same coefficients a as 
determined above. However, the extended matrix has one zero eigenvalue. The sum rule is 
just the inner product of the left null eigenvector xq with ^ = Xa, Xofi = x^Xa = 0. In 
particular, any set of moments which can be described by the initial seven-component chiral 
structure must satisfy the sum rule exactly. 

The sum rule is, in fact, satisfied rather accurately by the measured moments. The 
predicted value of /isa is /isa = (1-483 ± 0.012) /in, to be compared with the measured value 
I/WsaI = 1.610 ± 0.080) /iN- The discrepancy is at the level of 1.5cr or 8% of the measured 
transition moment. 

The situation is potentially different when loop corrections are included. While the ml^'^ 
corrections can be described using only the chiral structures considered so far and therefore 
satisfy the Okubo sum rule exactly [^], that is not the case for the corrections proportional 
to nislnms, a point misstated by those authors. The logarithmic corrections involve sums of 
terms of the form Mf ln(Mj/A) multiplied by coupling coefficients. If the mass dependence 
of the logarithms could be neglected, the sum rule would hold exactly for external meson 
masses that satisfy the Gell-Mann-Okubo mass formula. As a result, the chiral cutoff A 
drops out of the sum rule, and the violation depends only on the logarithm of the mass 
ratio, IuIMk/Mt^), hence would vanish for equal meson masses. However, the sum rule is 



fit to all eight measured quantities leads to essentially identical results because of the low weight 
accorded /isA in the fit. The unweighted fits used in [p|,|ll,10 distribute the uncertainty in ^sA 
over all the moments without changing the one test of the theory, as will be seen below. We believe 
that the present procedure is clearer. 



violated for Mk 7^ iVf,r, and the logarithmic loop corrections change the prediction for /i^a 
and introduce new structure in the 8x8 moment space.Q 

To investigate this point in more detail, we have constructed exact fits to the seven 
octet moments using the results for the nonanalytic parts of the loop corrections given in 
Eqs. (|3|-p!0D using the corrected coupling coefficients from 0, and have then calculated /isA- 
For illustration, we will present the results for the two cases considered in |^: (a) F = 0.4, 
D = 0.61, and C = —1.2, the one-loop axial vector couplings derived in |T^; and (b) F = 0.5, 



D = 0.75, and C = —1.5, the tree-level couplings. In both cases, we take the decup let-octet 
mass difference as 6 = 0.3 GeV, and use f^ = 0.093 GeV, fx = fr] = I-S/tt, /^t = —7.7, 
He = 1.94, and a chiral cutoff A = 1 GeV as in that reference. 

Although the sum rule is violated analytically, we find that the violation is remarkably 
small numerically. In particular, /i^a predicted to be 1.508 /xn for the parameters in case 
(a), rather the 1.483 /in as required by the sum rule. The difference is quite small compared 
to the experimental uncertainty in /isA- The results in case (b) are similar, with a predicted 
transition moment /i^a = 1-534. The chiral parameters for these two fits are given in Table 
II. 

It is evident from Table II that the fitted parameters fip and hd are consistent with 
SU(6): fio/f^F ~ 3/2. That is expected in the usual quark model picture. The axial 
coupling constants F, D, and C also satisfy the SU(6) relations. 

The result that loop corrections lead to only small violations of the Okubo sum rule 
appears to be quite robust. In particular, the changes in the loop corrections associated 
with different treatments of the meson-baryon couplings are largely absorbed by changes 
in the fitted counterterm parameters. We conclude that the theory with all the countert- 
erms present is not usefully predictive at the present time. The results do not, for example, 
distinguish in a useful way between the standard chiral approach with the nonanalytic con- 
tributions singled out, the approach of Meissner and Steininger |]lOl which includes analytic 



as well as nonanalytic terms, or the approach of Bos et al. [ITTI which allows only corrections 



that break SU(3) symmetry linearly, fits the data with just the seven parameters above, and 
satisfies the sum rule exactly. 

III. ANALYSIS OF THE LOOP CORRECTIONS 

The insensitivity of the overall fits to loop effects raises the question of whether there is 
any evidence that these effects are actually important. The difficulty is that the counterterm 
contributions are not calculable in the present theory, yet are essential for getting good 
overall fits to the data. A breakdown of the detailed contributions to the moments for case 
(b) is given in Table III. The results for case (a) are similar. 

The structure of the loop contributions is of considerable interest. We will examine this 
separately for the octet and decuplet cases. 



^The new structure can be parametrized by a single invariant Tr BMQBM of order 171"^ which 
affects only the A moment. 



A. Octet baryon contributions 

As shown in Table III, the contributions of type m}/'^ from the intermediate octet states 
are quite large, comparable to the tree-level contributions. The rrislnms corrections have 
the same general magnitude, but opposite signs. As a result there are large cancellations 
between these two types of nonanalytic contributions. This contradicts the arguments of 
chiral perturbation theory [0,|| and the l/Nc expansion 0J^J|,|] that the mglnms corrections 
should be small, of the same order of magnitude as the counterterms and ignorable relative 
to the my^ contributions. In fact, the smallness of the characteristic ratio |0 



V327r2/i: ^^ J \ Stt/J. / AnrriN A 

of the two types of nonanalytic terms, approximately -1/22 and -1/17 for X a pion or a kaon, 
is offset by very large ratios of the effective couplings, e.g., 16.1 and 12.7 for the pion and 
kaon contributions to the proton moment in case (b), and -28.6 and 9.9 for the S°, and the 
large cancellations ensue. There is no indication in this sector that the chiral perturbation 
series is converging. 

B. Decuplet contributions 

The contributions from intermediate decuplet states introduce a new complication, 
namely the dependence of the results on the octet-decuplet mass difference 6. The re- 
sults are very sensitive to 5, which is a QCD parameter rather than a chiral parameter.El 
Meissner and Steininger |T0[ argue that the decuplet contributions should be included in the 



chiral momentum expansion, expand in powers of q/S, and evaluate their couplings to take 
the resulting effects into account. Unfortunately, this argument requires that 6 be large on 
the scale of the momenta in the loop integrals. It is not, with M^^ < 6 < M^ for typical 
mass splittings 6 ~ 300 MeV, and the heavy decuplet argument fails. We conclude that 
the decuplet states must be regarded as light, on the same mass scale as the octet states, 
to obtain a consistent picture. Their contribution to the moments is then given by the ex- 
pressions in Eqs. (|3|-p!0D. The approximation in [^ becomes useful only for excited baryon 
multiplets, and can presumably be used to estimate those contributions. 

Jenkins et al. included the effects of the mass splitting in the m^^ terms, but took 5 = 
in the logarithmic contributions. This introduces very large changes in the latter, but does 
not affect their final conclusions because the differences are largely compensated through 
the counterterms. 

Using the value 5 = 300 MeV, we obtain the results for the decuplet contributions shown 
in Table III. The logarithmic contributions associated with the decuplet are large and tend 
to cancel the logarithmic corrections from the octet as suggested in [Q , but the cancellation 
is far from complete and all the logarithmic terms must be included. 



^E.g., 6 is described in the quark model in terms of a color hyperfine splitting between the spin-1/2 
octet and the spin-3/2 decuplet states. 



If we look at Table III from the point of view of a loop expansion, we see that the 
total contributions from the octet loops in Figs. 1 and 2 are considerably smaller than the 
individual m}-J'^ and mg In m^ terms because of the cancellations. This suggests that a loop 
expansion may converge more rapidly than the chiral expansion. However, the two types of 
terms tend to add rather than cancel for the decuplet, and turn out to dominate the final 
loop corrections. The counterterms are also quite significant. We see no convincing evidence 
that the chiral expansion with only the nonanalytic terms retained is under control. Our 
attempts to determine whether or not the loop corrections improve the tree-level fit to the 
moments are also inconclusive. We have refit the data omitting the counterterms, but with 
the loop corrections included and [Id and \ip readjusted. The quality of the fits is essentially 
unchanged with respect to tree level if all loop contributions are included. Jenkins et al. 
find that the fits can be improved by omitting the logarithmic corrections and readjusting 
F and D in addition to the somewhat small values quoted for case (a), but this procedure 
is not justified by the relative sizes of the terms retained and omitted (see Table III). 

C. Conclusions 

The foregoing results suggest that the loop expansion may give a more useful approach 
to the calculation of the baryon moments than the chiral and 1/A^c expansions. In particu- 
lar, the results obtained from the different diagrams in Figs. 1 and 2 are all of comparable 
size when combined with the relevant coupling factors and must be treated together to get 
meaningful results. However, loop effects are largely wiped out in the fits to the data when 
all the counterterms are treated as free parameters. Neither the chiral nor loop expansion 
is usefully predictive in this setting given the experimental and theoretical uncertainties. 
It is not possible, for example, to distinguish between the standard chiral expansion and 
the alternative approach proposed by Bos et al. |jTl|. We conclude that it will be neces- 



sary to develop a dynamical theory in which the baryons are treated as composite and the 
counterterms become calculable before the moment problem can be regarded as solved. 
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FIGURES 

1/2 

FIG. 1. The diagrams which give the non-analytic ms corrections to the baryon magnetic 
moments. Dashed Unes denotes mesons; single and double solid lines denote octet and decuplet 
baryons, respectively. 

FIG. 2. The diagrams which give the non-analytic ?Ti-sln?Ti-s corrections to the baryon magnetic 
moments. 
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TABLES 

TABLE L Table of the coefficients of the counterterm contributions to the baryon magnetic 
moments. This corresponds to the matrix X defined in the text. The notation is as fohows: 
oi = ^J'D^ 0,2 = fJ-Fj Eq- (1)) while 03, . . . , 07 are the parameters are defined in Eq. (2). The S" 
moment is not included, but is given by the SU(2) relation ^^o = (/"s+ + /^E-)/2- 

Baryon oi 02 03 04 05 uq a-j 



1 1 1 1 1 1-1 

3 3 3 3 

-- -- -- — - 

3 " " " 3 3 3 

i 1 -i 

i -1 -i 



1 


1 


1 











1 








-1 

















-1 


1 


-1 





















V 
n 

s- 

-° — - - — - — i 

" 3UUU3 33 

-- i -1 1 -1 -i i -i 

" 3^-^^33 3 

A -i -| -| 

TPK ^0000 



1 


1 


3 


3 


2 


2 


3 


3 














2 


2 


3 


3 


1 


1 


3 


3 





8 
~9 









TABLE IL The best fit parameters obtained by fitting the seven measured octet magnetic 
moments exactly. The parameters are given in units of ^n- The SA transition moment is predicted 
to be ;UsA = 1.508^N in case (a), 1.534^n in case (b), and (1.483ib0.012) ^un by the Okubo relation, 
case (c). 



Case 


fli 


^2 


as 


a4 


05 


ag 


a? 


(a) 


3.183 


L992 


0.269 


-0.215 


0.534 


0.224 


-0.392 


(b) 


3.946 


2.353 


-0.001 


-0.172 


0.569 


0.694 


-1.165 


(c) 


2.568 


L809 


0.591 


-0.253 


0.497 


-0.506 


0.621 



TABLE in. Detailed breakdown of the contributions of the loop integrals and the counterterms 
to the fitted magnetic moments of the octet baryons for the parameters of case (b). Table IL The 
first seven moments were used as input, and are fitted exactly. The only prediction of the theory 
is /USA- 



MS 


I^D, I^F 


m'J\N) 


mV\A) 


In ms{N) 


Inms(A) 


Loops 


CT 


MB 


p 


3.668 


-2.227 


-0.162 


2.148 


-1.271 


-1.512 


0.636 


2.793 


n 


-2.631 


0.832 


0.573 


-1.491 


1.258 


1.172 


-0.454 


-1.913 


S+ 


3.668 


-2.843 


-0.474 


2.829 


-1.111 


-1.598 


0.388 


2.458 


s- 


-1.037 


0.611 


-0.349 


-0.804 


0.031 


-0.511 


0.388 


-1.160 


so 


1.315 


-1.116 


-0.411 


1.012 


-0.540 


-1.055 


0.388 


0.649 


A 


-1.315 


1.116 


0.411 


-1.211 


0.614 


0.931 


-0.228 


-0.613 


rP 


-2.631 


2.288 


0.698 


-2.650 


0.739 


1.075 


0.306 


-1.250 


L^ 


-1.037 


1.339 


-0.286 


-1.324 


0.057 


-0.215 


0.601 


-0.651 


yPk 


2.278 


-1.157 


-0.496 


1.680 


-0.771 


-0.745 


0.000 


1.534 
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Fig. 1 



(b) 
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Fig. 2 



